Langevin dynamics of the pure SU{2) deconfining transition 



E. S. FragaI*, G. KREiN^t, A. J. MizherI* 

^Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68528, 21941-972 Rio de Janeiro, RJ , Brazil 
Instituto de Ftsica Teorica, Universidade Estadual Paulista, 
Rua Pamplona 145, Sao Paulo, SP 01405-900, Brazil 



o 
o 

(N 

< 
O 

m 



Oh 
I 

Oh 

CD 



(N 
> 

(N 
p 

o 
o 



X 



We investigate the dissipative real-time evolution of the order parameter for the deconfining tran- 
sition in the pure SU{2) gauge theory. The approach to equilibrium after a quench to temperatures 
well above the critical one is described by a Langevin equation. To fix completely the markovian 
Langevin dynamics we choose the dissipation coefficient, that is a function of the temperature, 
guided by preliminary Monte Carlo simulations for various temperatures. Assuming a relationship 
between Monte Carlo time and real time, we estimate the delay in thermalization brought about by 
dissipation and noise. 



I. INTRODUCTION 

The study of the dynamics of phase conversion dur- 
ing the deconfinement transition for pure gauge theories 
might shed some light into the process of thermalization 
of the quark-gluon plasma in hot QCD in a controlled 
fashion. Indeed, for pure SU{N), that can be seen as 
QCD in the limit of infinitely heavy quarks, lattice sim- 
ulations are well developed, yielding a precise prediction 
for the deconfinement critical temperature and a good 
understanding of the corresponding thermodynamics 1] . 
In this limit there is a global Z{N) symmetry associated 
with the center of the gauge group, so that one can use 
the Polyakov loop to construct a well-defined exact order 
parameter 0, 0] , and an effective Landau-Ginzburg field 
theory based on this quantity 0, @ . 

The effective potential for T << Td, where Td is the 
deconfinement critical temperature, has only one mini- 
mum, at the origin, where the whole system is localized. 
With the increase of the temperature, N new minima ap- 
pear. At the critical temperature, T^, all the minima are 
degenerate, and above the new minima become the 
true vacuum states of the theory, so that the extremum 
at zero becomes unstable or metastable and the system 
starts to decay. In the case of SU{3), whithin a range of 
temperatures close to there is a small barrier due to 
the weak first-order nature of the transition [gJ , and the 
process of phase conversion will thus be guided by bub- 
ble nucleation. For larger T, the barrier disappears and 
the system explodes in the process of spinodal decompo- 
sition Q. For SU{2), the transition is second-order]^, 
and there is never a barrier to overcome. 

Real-time relaxation to equilibrium after a thermal 
quench followed by a phase transition, as considered 
above, can in general be described by standard reaction- 
diffusion equations For a non-conserved order pa- 
rameter, ip{x,t), such as in the case of the deconfining 
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transition in pure gauge theories, the evolution is given 
by the Langevin, or time-dependent Landau-Ginzburg, 
equation 



dip _ SF 
'dt ~ ^Jii 
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dip 



(1) 



where F = F{ip,T) is the coarse-grained free energy of 
the system, a is the surface tension and U = U{ip,T) is 
the effective potential. The quantity F is known as the 
dissipation coefficient and will play an important role in 
our discussion. Its inverse defines a time scale for the 
system, and is usually taken to be either constant or as 
a function of temperature only, F — F(T). The function 
^ is a stochastic noise assumed to be gaussian and white, 
so that 
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according to the fluctuation-dissipation theorem. From 
a microscopic point of view, the noise and dissipation 
terms are originated from thermal and quantum fluctu- 
ations resulting either from self-interactions of the field 
representing the order parameter or from the coupling 
of the order parameter to different fields in the system. 
In general, though, Langevin equations derived from a 
microscopic field theory [9;] contain also the influence of 
multiplicative noise and memory kernels [l^, [ill, • 

In this paper, we consider the pure SU (2) gauge the- 
ory, without dynamical quarks, that is rapidly driven to 
very high temperatures, well above T^, and decays to the 
deconfincd phase via spinodal decomposition. We are 
particularly interested in the effect of noise and dissipa- 
tion on the time scales involved in this "decay process" , 
since this might provide some insight into the problem 
of thermalization of the quark-gluon plasma presumably 
formed in high-energy heavy ion collisions [13|. For the 
order parameter and effective potential we adopt the ef- 
fective model proposed in Ref. and the choice of the 
dissipation coefficient, that is a function of the temper- 
ature, is guided by preliminary Monte Carlo simulations 
for various temperatures, comparing the short-time ex- 
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ponential growth of the two-point Polyakov loop corre- 
lation function predicted by the simulations to the 
Langevin description assuming, of course, that both dy- 
namics are the same (see, also, the extensive studies of 
Glauber evolution by Berg et al. [IBl)- This procedure 
fixes completely the Markovian Langevin dynamics, as 
will be described below, if one assumes a relationship be- 
tween Monte Carlo time and real time. Once the setup 
is defined for the real-time evolution, we can estimate 
the delay in thermalization brought about by dissipation 
and noise by performing numerical calculations for the 
dynamics of the order parameter on a cubic lattice. As 
will be shown in the following, the effects of dissipation 
and noise significantly delay the thermalization process 
for any physical choice of the parameters, a result that 
is in line with but is even more remarkable than the one 
found for the chiral transition p^ . 

The paper is organized as follows. In Section II, we de- 
scribe the effective model adopted for the Langevin evo- 
lution implementation, as well as the analytic behavior 
for early times. In Section III, we discuss the necessity of 
performing a lattice renormalization to have results that 
thermalize to values that are independent of the lattice 
spacing and free from ultraviolet divergences, and present 
the necessary counterterms. In Section IV we briefly de- 
scribe the Glauber dynamics of pure lattice gauge theory 
that can be used to extract the dissipation coefficient for 
different values of the temperature. Details and quanti- 
tative results from lattice simulations will be presented 
in a future publication [17J . In Section V we present and 
analyze our numerical results for the time evolution of 
the order parameter for deconfinement after a quench to 
temperatures above T^. Finally, Section VI contains our 
conclusions and outlook. 



II. EFFECTIVE MODEL AND LANGEVIN 
DYNAMICS 

Since we focus our investigation on pure gauge SU (N) 
theories, we can adopt effective models built by using 
functions of the Polyakov loop as the order parameter for 
the deconfining phase transition. If quarks were included 
in the theory, the Z{N) symmetry present in pure glue 
systems would be explicitly broken, and the Polyakov 
loop would provide only an approximate order parame- 
ter. For euclidean gauge theories at finite temperature, 
one defines the Polyakov loop as: 



P{x) = r exp 



«5 



l/T 



drAo {x, t) 



(3) 



where T stands for euclidean time ordering, g is the gauge 
coupling constant and Aq is the time component of the 
vector potential. 

The effective theory we adopt [3] is based on a mean- 
field treatment in which the Polyakov loops are constant 
throughout the space. The degrees of freedom that will 



be used to construct the free energy are the eigenvalues 
of the Polyakov loop, rather than {TrpP{x)). Working 
in SU (N) gauge theories the Polyakov loop is unitary, so 
that it can be diagonalized by a unitary transformation, 
assuming the form 



Pjk = exp(i6lj) S-jk 



(4) 



At one loop, the free energy for gluons in a constant 
Aq background is given by 



fpertiO) - ln[det{-DlaA] , 



(5) 



where Dadj is the covariant derivative acting on fields 
in the adjoint representation. This expression can be 
written in a more explicit form: 



(6) 



where B is defined in Eq.(|4]), and A^j^ = 9j — 9k- Here 
we have the "bare" dispersion relation = |k|. In or- 
der to include confinement in this effective model de- 
scription, one can introduce an ad hoc "thermal mass" 
for the gluons, so that the dispersion relation becomes 
cjk = vk? + M^ . The value of M can be related to 
the critical temperature Td extracted from lattice simu- 
lations. 

Parametrizing the diagonalized Polyakov loop as 
dia(7[exp(i0^/2j ••■J ^'/'ij "^i^ii ■■■7 ~*</'Af/2]j we can con- 
struct the effective potential from the free energy above. 
For SU{2), it can be written in the following convenient 
form: 



(7) 



where we have defined -0 = 1 — 20/7r, and used the re- 
lation between the mass M and the critical temperature 
Td — a/SM^T^tt^- In Fig- Uwe display U as a function 
of V' for different values of the temperature. One can 
see from this plot that for T <C the minimum is at 
'(/'(x, t) = 0. As the temperature increases new minima 
appear; above the critical temperature they become the 
true vacuum states of the system. Now, if at t = the 
temperature is rapidly increased to T ^ Td, the system 
is brought to an unstable state and therefore will start 
"rolling down" to the new minima of the effective poten- 
tial. 

To study the time evolution, we consider a system char- 
acterized by a coarse-grained free energy of the form 



F{'iP,T)= d^x -(VV-)^ + C/(V',T) 



(8) 



where U is the effective potential obtained above, and a 
plays the role of a surface tension [l^, assuming the fol- 
lowing value for SU{2): a = TT^T/g^, g being the gauge 
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FIG. 1: Effective potential for tfie case of SU{2) for different 
values of the temperature. 



coupling. The approach to equilibrium of the order pa- 
rameter ipi^j i) will then be dictated by the Langevin 
equation ([1]) that, for arbitrary times, has to be solved 
numerically on a lattice. 

At very short times, however, when ^(x, t) « 0, non- 
linear terms in the evolution equation can be neglected, 
so that Eq. ^ reduces to 



(9) 



in the Fourier space, where mx is a dimensionless thermal 
mass that can be written as 



(10) 



6 is a number that depends on the details of the quadratic 
term of the particular effective potential adopted, so that 
it will be different, for instance, if we consider SU{2) 
{b = 1) or SU{3) {b — 10/9). One can, then, approximate 
the (noiseless) solution in Fourier space by ip(k, t w 0) ~ 
e^^* , where Ak are the roots of the quadratic equation 



A£ + ( - 1 A, 



k2- 



= 0. 



For wavelength modes such that 
|k| < 



\ a 



1/2 



(11) 



(12) 



one has already oscillations, but those are still damped 
by a factor exp(— ri/2CT). It is only for longer wavelength 
modes, i.e. 



|k| <k,^ 



(13) 



that there will be an explosive exponential growth corre- 
sponding to the regime of spinodal decomposition. 



As time increases, the order parameter increases and 
non-linear contributions take over. To study the com- 
plete evolution of the phase conversion process, we have 
to solve (m numerically on a lattice. In the next section 
we discuss the need for lattice renormalization to avoid 
spurious ultraviolet divergences in the dynamics. 



III. LATTICE RENORMALIZATION 

In performing lattice simulations of the Langevin evo- 
lution, one should be careful in preserving the lattice-size 
independence of the results, especially when one is con- 
cerned about the behavior of the system in the continuum 
limit. In fact, in the presence of thermal noise, short and 
long wavelength modes are mixed during the dynamics, 
yielding an unphysical lattice size sensitivity. The issue 
of obtaining robust results, as well as the correct ultra- 
violet behavior, in performin g L angevin dynamics was 
discussed by several authors [IQI . uuL [21| . [22 . [23j . The 
problem, which is not a priori evident in the Langevin 
formulation, is related to the well-known Rayleigh- Jeans 
ultraviolet catastrophe in classical field theory. The dy- 
namics dictated by Eq. ([1]) is classical, and is ill-defined 
for very large momenta. 

Equilibrium solutions of the Langevin equation that 
are insensitive to lattice spacing can be obtained, in prac- 
tice, by adding finite-temperature counterterms to the 
original effective potential, which guarantees the correct 
short-wavelength behavior of the discrete theory. Fur- 
thermore, it assures that the system will evolve to the 
correct quantum state at large times. For a more de- 
tailed analysis of lattice renormalization in the Langevin 
evolution, including the case of multiplicative noise, see 
Ref. [241 . 

Since the classical scalar theory in three spatial di- 
mensions is super-renormalizable, only two Feynman di- 
agrams are divergent, the tadpole and sunset. The sin- 
gular part of these graphs can be isolated using lattice 
regularization, and then subtracted from the effective po- 
tential in the Langevin equation. For a scalar field theory, 
explicit expressions for the counterterms were obtained 
by Farakos et al. [2^ within the framework of dimen- 
sional reduction in a different context. 

Following Ref. [1^, we write the bare potential in a 
three-dimensional field theory in the following form 

V(0) = -im^^^ ^ iA30^ (14) 



where m is the bare mass of the field (f> and the subindex 
in A3 stresses the fact that this is the coupling of a three- 
dimensional theory. In Ref. [1^, this dimensionally- 
reduced theory was obtained from a four-dimensional 
theory with a dimensionless coupling A, assuming a 
regime of very high temperature. At leading order, one 
has A3 = A T. The mass counterterm, which is defined 
such that 



1 



2 i2 

m (j) 
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(15) 
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is given by 



6m^ = 3 X: 



0.252731 



6X1 



ln( A 



0.09 



(16) 

where a is the lattice spacing and fi is the renormaUzation 
scale. The first term conies from the tadpode diagram 
and the second one from the sunset. Finite constants are 
obtained imposing that, after renormahzation, the sunset 
diagram yields the same value for three renormahzation 
schemes: lattice, momentum subtraction and MS [2^ . 
Notice that in order to obtain lattice-inde pen dent results 
physical quantities become /u-dependent [2l|. However, 
since the contribution from the /i-dependent term is log- 
arithmic, variations around a given choice for this scale 
affect the final results by a numerically negligible factor, 
as we verified in our simulations, so that this dependence 
is very mild. 

Since the field tp in the effective model we consider 
here is dimensionsless, it is convenient to define the di- 
mensionful field f = a^/'^ip in order to relate results from 
Ref. [25| to our case more directly. 

Now we can write our Langevin equation, Eq. Ilj, in 
terms of the field tp. For SU{2), we have 



^9lP 2 , X 3 ^ 



where 



the theory is super-renormalizable and only mass coun- 
terterms are needed. In terms of the original ip, our renor- 
malized Langevin equation is finally given by 



r 



'dt 
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where 
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(23) 
(24) 



One can factor out the appropriate powers of T in this 
expression to make explicit the mass dimensions: 



5Ml 



3T- 



6r^ 



0.252731 (g 



167r2 
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0.09 



(25) 



Notice that for sufficiently high temperatures the sym- 
metry of the potential is restored. 



IV. DISSIPATION COEFFICIENT FROM 
MONTE CARLO EVOLUTION 



Xl = 



T^ml 
a 

3<t2 



e(T) 



(18) 
(19) 



The subindex L in these quantities is a reminder that 
they refer to the Langevin equation. It is clear that 
Eq. (fT7| corresponds to an effective action S{ip) given 
by 



1 



Siip) = -(V^) 



1 



2 2 



1 



(20) 



Once we have identified the mass term and the cou- 
pling constant, we can renormalize the Langevin equa- 
tion, which becomes 



Tdp 



Xl'P' 



2 ' Sml)tp 



(21) 



where 



3Xr 



0.252731 



6Af 



IGtt^ 







['"(-) 


+ 0.09 


\^J.aJ 





(22) 

Notice that we have used the same symbol (p to denote 
both the renormalized and non-renormalized fields, since 



In lattice simulations for pure SU (N) gauge theories, 
one can implement the Glauber dynamics by starting 
from thermalized gauge field configurations at a temper- 
ature T < Td and then changing the temperature of the 
entire lattice that is quenched to T > [l^, [l3| . The 
gauge fields are then updated using the heat-bath algo- 
rithm of Ref. [1^ without over-relaxation. A "time" unit 
in this evolution is defined as one update of the entire 
lattice by visiting systematically one site at a time. 

The structure function, defined as 



S{k,T) = {LF{k,T)LFhk,T)) 



(26) 



where Lp{k,T) is the Fourier transform of Lf(x,t), the 
Polyakov loop in the fundamental representation, can be 
used to obtain the values of the dissipation coefficient, F, 
for different values of the final temperature, T, as follows. 
At early times, immediately after the quench, ~ and 
one can neglect the terms proportional to and "0^ in 
the effective potential to first approximation. It is not 
difficult to show that at early times, when ^ is small, the 
structure function can be written as 



S{k,T) = S{k,Q) exp[2cj(fc)T] , 



where 



uj{k) 



ikl - 



(27) 



(28) 
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In obtaining this expression we have neglected the 
second-order time derivative in Eq. Q, which should 
be a good approximation for a rough estimate of T. For 
the effective potential adopted here, is given by 



we compute 



9M^ 
1^ 



(29) 



One sees that for momenta smaller than the critical mo- 
mentum fcc, one has the familiar exponential growth, sig- 
naling spinodal decomposition. Plotting lnS{k,T)/T for 
different values of k allows one to extract 2uj{k) and, in 
particular, the value of k^. Once one has extracted these 
values, r can be obtained from the following relation: 



= uj{k) 



9 



7r2r(fc2 _ fc2) • 



(30) 



Now, in Monte Carlo simulations one does not have a 
time variable in physical units and so, by plotting In S 
from the lattice, one obtains values of 2uj{k) that do not 
include the (unknown) scale connecting real time r and 
Monte Carlo time. Nevertheless, if one assumes that the 
relation between the Langevin time variable r and the 
Monte Carlo time is linear, one can parametrize this re- 
lation in terms of the lattice spacing a as r = a\MCi 
where \mc is a dimensionless parameter that gives this 
relation in units of the lattice spacing. An estimate for 
the relationship between Monte Carlo time and real time 
is given in Ref. 27], where the authors evaluate the num- 
ber of sweeps necessary for the system to freeze-out. In 
this reference, the authors implement lattice Monte Carlo 
simulations of the change of the Polyakov loop under lat- 
tice expansion and temperature falloff. The freeze-out 
number of sweeps was defined as being the number of 
sweeps necessary for the Polyakov loop to reach zero. 
This number was found to be of the order of 5000 for the 
range of temperatures we are considering here. Using the 
phenomenological value of 9 fm/c [l^ as the freeze-out 
time, one can then obtain \mc- 

Preliminary simulations clearly show that F"^ de- 
creases as the final temperature increases Guided 
by these results, we choose in the case of SU{2) T{T) — 
10'^ fm~2 for our Langevin simulations, which we describe 
in the next section. 



V. NUMERICAL RESULTS FOR 
DECONFINEMENT AND DISCUSSION 

We solve Eq. ([1]) numericaly for SU{2) in a cubic 
spacelike lattice with 64'^ sites under periodic bound- 
ary conditions, using the semi-implicit finite-difference 
method for time discretization and finite-difference Fast 
Fourier Transform for spatial discretization and evolu- 
tion [2^ . To compute the expectation value of the order 
parameter ip^ we average over thousands of realizations 
with different initial conditions around tp and differ- 
ent initial configurations for the noise. At each time step 



it), 



(31) 



ijk 



where the indices i,j,k indicate the position of the site 
on the lattice. 

The thermal mass M can be determined through the 
deconfinement temperature. For SU{2), Td — 302 MeV, 
so that M = 775 MeV. In Fig. [2] we show the time 
evolution of {tp) for the SU{2) case, normalized by i/'Oj 
which corresponds to the value of the order parameter 
at the vacuum. The dotted line represents the case with 
no noise and no dissipation, the dashed line corresponds 
to the case with only dissipation, and the full line to the 
complete case. Simulations were run under a tempera- 
ture of r — 6.6 Td, which ensures that there is no barrier 
to overcome, and the dynamics will be that of spinodal 
decomposition. For this temperature the value of F is 
given by 10'' fm^^^ jj-i accordance with the discussion of 
the previous section. 



1.5 
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FIG. 2: Langevin evolution of the SU{2) order parameter. 



One can clearly see from the figure that dissipation 
brings strong effects in the time evolution, delaying con- 
siderably the necessary time for the onset of the decay 
process. Noise acts in the same direction as dissipa- 
tion, retarding even more the time of equilibration: from 
around 2 fm/c, for the simulation including dissipation 
effects only, to more than 4 fm/c in the complete case. 
Comparing our results to those from a similar calculation 
performed for the case of the chiral phase transition [16| , 
it is evident that in the former dissipation and noise have 
similar but stronger effects. This might signal that the 
dynamics of the deconfinement transition is more sensi- 
tive to medium effects. However, this is a very prema- 
ture conjecture, since both effective theory approaches 
are rather simplified descriptions of in- medium QCD. 
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VI. CONCLUSIONS AND OUTLOOK 

We have presented a systematic procedure to study 
the real-time dynamics of pure gauge deconfinement 
phase transitions, considering in detail the case of SU{2). 
Given an effective field theory for the order parameter 
of the transition, we have discussed the necessity to in- 
troduce counterterms from lattice renormalization that 
guarantee lattice independence of physical results. These 
counterterms were computed for the case of SU{2) or 
any theory whose effective model exhibits the same di- 
vergence structure. 

For the Langevin evolution, one needs the dissipation 
coefficient as an input. We have described a recipe to 
extract this kinetic quantity from Glauber dynamics in 
Monte Carlo simulations. The value adopted here is 
based on preliminary lattice results. A detailed analysis 
will be presented in a future publication (l7j . together 
with Langevin evolution results for the case of SU{3). 

From our results for the dynamics of the deconfining 
transition in SU{2), we conclude that dissipation and 
noise play a very relevant role, being responsible for de- 
lays in the equilibration time of the order of 100%. So, 
effects from the medium are clearly significant in the de- 
termination of the physical time scales, and should be 
included in any description. 

Of course, the treatment implemented here is very sim- 
plified in many respects. First, there is a need for a more 
robust effective theory for the order parameter of the de- 
confining transition. Recently, studies of the renormal- 
ization of Polyakov loops naturally lead to effective ma- 
trix models for the deconfinement transition [SO^ , unfold- 



ing a much richer set of possibilities than the approach 
considered here. In particular, eigenvalue repulsion from 
the Vandermonde determinant in the measure seems to 
play a key role as discussed in Ref. [3l|. Nevertheless, 
these studies have shown that, in the neighborhood of 
the transition, the relevant quantity is still the trace of 
the Polyakov loop. 

Second, there is a need to construct a phenomeno- 
logical generalized Landau-Ginzburg effective theory de- 
scribing simultaneously the processes of chiral symmetry 
restoration and deconfinement in the presence of massive 
quarks as discussed in Ref. [s^l- Then, the dynamics of 
the approximate order parameters, the chiral condensate 
and the expectation value of the trace of the Polyakov 
loop, will be entangled. Finally, if one has the physics 
of heavy ion collisions in mind, effects brought about by 
the expansion of the plasma [33] and by its finite size |34| 
will also bring corrections to this picture. 

In a more realistic approach, time scales extracted from 
the real-time evolution of the order parameters can be 
confronted with high-energy heavy ion collisions exper- 
imental data, and perhaps provide some clues for the 
understanding of the mechanism of equilibration of the 
quark-gluon plasma presumably formed at Relativistic 
Heavy Ion Collider (RHIC). 
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